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Plasma turbulence is studied using Doppler reflectometry at the TJ-II stellarator. By scanning 
the tilt angle of the probing beam, different values of the perpendicular wave numbers are probed 
at the reflection layer. In this way, the interaction between zonal flows and turbulence is reported 
with (a) spatial, (b) temporal, and (c) wavenumber resolution for the first time in any magnetic 
confinement fusion device. 

We report measurements of the bicoherence across the Low to High (L-H) confinement transition 
at TJ-II. We examine both fast transitions and slow transitions characterized by an intermediate (I) 
phase. The bicoherence, understood to reflect the non-linear coupling between the perpendicular 
velocity (zonal flow) and turbulence amplitude, is significantly enhanced in a time window of several 
tens of ms around the time of the L-H transition. It is found to peak at a specific radial position 
(slightly inward from the radial electric field shear layer in H mode), and is associated with a 
specific perpendicular wave number (k± ^6 — 12 cm -1 ). In all cases, the bicoherence is due to the 
interaction between high frequencies and a rather low frequency, as expected for a zonal flow. 



The spontaneous transition to the High confinement 
regime in fusion plasmas is a topic of great interest, due 
to its potential advantages for the design of a fusion reac- 
tor and the complex nonlinear physics involved. At the 
TJ-II stellarator, the transition from Low to High con- 
finement (L-H) is "soft" in the sense that the confine- 
ment improvement factor is small, and yet it possesses 
the typical features of any L-H transition (rapid drop of 
H a emission, reduction of fluctuation amplitudes, forma- 
tion of a sheared flow layer, ELM-like bursts) [HQ. 

A common explanation for the L-H transition involves 
the generation of a zonal flow in a narrow radial region 
of the plasma near its edge. The zonal flow would be 
driven by turbulence via the Reynolds Stress mechanism 
in a strong pressure gradient region. The growth of the 
(fluctuating) zonal flow then leads to the suppression of 
turbulence, creating a transport barrier, and the subse- 
quent establishment of a steady state sheared flow with 
associated radial electric field that maintains the bar- 
rier [3|— . The detailed observation of this sequence of 
events is a challenge. 

Due to the temporal and spatial scales involved in the 
L-H transition physics, specific experimental techniques 
are required to investigate the turbulence and flow dy- 
namics 0. At TJ-II, a two-channel Doppler reflectome- 
ter is available, allowing the simultaneous measurement 
of the perpendicular plasma velocity and density fluctua- 
tions at two radial positions, with good spatial and tem- 
poral resolution. In addition, the Doppler reflectometer 
offers the possibility to probe different turbulence scales 
(wave vectors) by steering the probing beam 0, Q . 

Experimental setup - The TJ-II vacuum magnetic 
geometry is completely determined by the currents flow- 
ing in four external coil sets. At TJ-II, the normalized 
pressure (/?) is generally low, even in discharges with Neu- 



tral Beam Injection (NBI), and currents flowing inside 
the plasma are generally quite small (unless explicitly 
driven), so that the actual magnetic configuration and 
the rotational transform l/2tt are typically rather close 
to the vacuum magnetic configuration [8] . Here we focus 
on two magnetic configurations having i(a)/27T = 1.630 
and 1.553. 

The experiments have been carried out in pure NBI 
heated plasmas (line averaged plasma density (n e ) = 2 — 
4x 10 19 m -3 , central electron temperature T e = 300 — 400 
eV). The NBI input heating power is kept constant at 
about 500 kW during the discharge but the fraction of 
NBI absorbed power - taking into account shine through, 
CX and ion losses, as estimated using a simulation code 
- increases from 55 to 70% as the plasma density rises. 

In this work, we consider two different L-H transition 
scenarios. Standard or fast L-H transitions are observed 
in the magnetic configuration having i(a)/2n = 1.630 
while slow transitions are observed in the configura- 
tion having l(o)/2k = 1.553 and a low order rational 
(3/2) at p = r/a ~ 0.72. In the latter configuration, 
the so-called intermediate phase (I), characterized by a 
predator-prey type interaction between turbulence and 
flows, appears between the L and H phases (9MT1 ] . which 
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and on other de- 



has also been seen in models 

vices [3, [Hj. In both scenarios, spatiotemporal and scale 
resolved Doppler reflectometry measurements were per- 
formed in series of repetitive discharges 0, 0] • 

Doppler reflectometry - In Doppler reflectometry, 
a finite tilt angle is purposely introduced between the 
incident probing beam and the normal to the reflecting 
cut-off layer, and the Bragg back-scattered signal is mea- 
sured [16|. The amplitude of the recorded signal, A, is 
a measure of the intensity of the density fluctuations, ft. 
By scanning the tilt angle of the probing beam, different 
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perpendicular wave numbers are probed, thus allowing 
scale-resolved turbulence measurements. Furthermore, 
as the plasma rotates in the reflecting plane (flux sur- 
face), the scattered signal experiences a Doppler shift. 
The size of this shift is directly proportional to the rota- 
tion velocity of the plasma turbulence perpendicular to 
the magnetic field lines, v±, and therefore to the plasma 
background ExB velocity, provided the latter dominates 
over the phase velocity of density fluctuations (cf. [2[). 
The Doppler reflectometer signals, sampled at 10 MHz, 
allow determining h and v± with high temporal and spa- 
tial resolution. 

The relation between zonal flow growth, involving an 
interaction between high frequencies associated with drift 
wave turbulence and low frequencies associated with the 
zonal flow, and the cross bispectrum was elucidated 
in [I?}. Probing zonal flow dynamics therefore requires 
measuring the cross bispectrum between the radial and 
poloidal ExB velocities of high-frequency drift waves 
and the low-frequency zonal flow. An early attempt at 
observing the bicoherence across the L-H transition was 
made at DIII-D [lij . As such measurements are quite dif- 
ficult to perform in the plasma interior, the bicoherence 
is often calculated using density fluctuations only, which 
complicates the interpretation of the results. However, 
Doppler reflectometry allows the measurement of the tur- 
bulence fluctuation amplitude h as well as the fluctuating 
perpendicular flow v± , both with good temporal and spa- 
tial resolution, making the two main quantities involved 
in zonal flow dynamics accessible experimentally. In this 
work, we consider the complex amplitude Doppler reflec- 
tometry signal, Ae 1 ^ = (A cos 0, Asmc/)), which contains 
information of both the perpendicular flow, v± = d<j)/dt, 
and the density turbulence, n oc A. 

Bicoherence analysis - The cross bicoherence of two 
signals f(t) and g(t) is defined as [19| 
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where f(uS) is the Fourier transform of /(£), a star in- 
dicates the complex conjugate, and the angular brackets 
refer to an average over equivalent realizations. In prac- 
tice, the equivalent realizations are successive time sec- 
tions, assuming the system is in steady state. The value 
of b 2 is bounded between and 1, so that a high value 
indicates wave-wave coupling between the two base fre- 
quencies uji and Co>2 and their sum frequency, uo\ + GJ2- 
The auto-bicoherence is b 2 



We also report the "summed bicoherence", defined as 
frsiimM = (& S ^1^2)) Wl + W2 =u;, and the "total bicoher- 
ence", defined as 6 2 otal = {b 2 (uj\, uj 2)) oj u oj 2 ' To judge the 
relevance of a bicoherence value, it is compared to the ex- 
pectation value of the bicoherence under the assumption 
that the signal is uncorrelated random noise. 

Computing the auto-bicoherence for the complex 
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FIG. 1: Discharge 22277. Configuration with i{a)/2n = 
1.704. Top: Mean auto-bicoherence and mean summed auto- 
bicoherence of complex Doppler reflectometry channel 1. Bot- 
tom: Cross bicoherence between v± and h (Doppler reflectom- 
etry channel 1) during the L-H transition (1100 < t < 1130 
ms, tL-H = 1116 ms). The dashed line indicates the noise 
level. 



Doppler reflectometer data around the time of the L- 
H transition, one typically obtains results such as those 
shown in Fig.[TJ The auto-bicoherence is very significant, 
as the value is some two orders of magnitude above the 
noise level. The auto-bicoherence is concentrated along 
the line fa — —fi/2 (the 'first harmonics' line), fa — fi 
(interaction with a very small difference frequency), and 
also at fa — 0. These are indications that high frequen- 
cies / a , fa of the order of 1-2 MHz are interacting via a 
very small difference frequency f c — fa~ fb (or harmon- 
ics fa = 2f a - fb) such that \f c \ < \f a \, \f b \. Frequencies 
above 1 MHz are not commonly observed in this context 
due to technical limitations; however, note that similar 
results were obtained at DIII-D [18|. 

To facilitate the interpretation of the auto-bicoherence, 
Fig. [1] also shows the cross bicoherence between v±_ and 
n for the same discharge. Here, significant coupling is 
limited mainly to the lines fi — and fa c± indicat- 
ing an interaction between two very similar frequencies 
f a ,fb, interacting via a very small difference frequency 
fc = fa~ fb such that \f c \ <C |/ |, \f b \, while harmonics 
are weak. This is in fact what one would expect for the 
interaction between a zonal flow and turbulence 
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The small frequency of the triplet is at or below the fre- 
quency resolution (50 kHz), whereas the two high fre- 
quencies are in the range from 1 to 2 MHz. The fact 
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that the cross bicoherence graph shows linear structures, 
while the auto-bicoherence exhibits additional harmonic 
structures, is a clear indication that the physical inter- 
action primarily involves v± and h. Although the cross 
bicoherence between v± and h is easier to interpret, the 
auto-bicoherence of the complex Doppler reflectometry 
signal (containing the same information) is easier to com- 
pute, not affected by the numerical errors inherent in the 
velocity calculation, and yields a larger signal to noise ra- 
tio in the total bicoherence value. Therefore, we will use 
the latter as a generic proxy to study the spatiotemporal 
evolution of the total bicoherence. 

The fast L— H transition - To study the radial and 
temporal behavior of the bicoherence across the fast L-H 
transition, we have analyzed a series of 25 discharges in 
the same magnetic configuration (L(a)/2ir = 1.630). See 
Ref. [6] for more details about this series of discharges. 
The Doppler reflectometer probing frequency and tilt an- 
gle were varied on a shot to shot basis, and correspond- 
ingly the radial position of the reflecting layer and the 
wave number k±. Fig. [2^ indicates the evolution of the 
density profile across the L-H transition. Due to the 
evolving profile, both the position of the reflecting layer 
and the wavenumber change as a function of time across 
the transition jgj. The measured channel locations and 
wave numbers are shown in Fig. [2Jd. Note that k± varies 
by less than 5% between the L and H phases (matching 
black and red dots). The p and k± values indicated on 
the axes of Figs. [2t-e correspond to the situation in H- 
mode. Note that p is a vacuum flux surface label; the ac- 
tual confinement radius is somewhat less than p — 1 [211. 



The temporal resolution of the bicoherence calculation is 
2 ms. 

Figs. [2t-e show the values of the auto-bicoherence 
b 2 (k±,p, At), averaged over k±, p, and At, respectively. 
It is seen that the bicoherence appears ~ 20 — 30 ms be- 
fore the L-H transition at a specific radius (p < 0.8) and 
that specific perpendicular wave numbers are involved 
(fc_i_ ~ 8 — 11 cm -1 ). The bicoherence persists until about 
10 — 20 ms after the L-H transition. Similar behavior is 
observed across a wide range of magnetic configurations, 
showing that the phenomena are not associated with any 
specific rational surface. 

It should be noted that turbulence is reduced at all 
scales, but most in intermediate scales, k± ~ 7—11 cm -1 , 
as reported in earlier work [6]. After the transition, the 
radial electric field shear layer is located at p ~ 0.83, 
somewhat outward from the position of the maximum 
observed bicoherence values. 

The slow L— I— H transition - Fig. [3] shows the auto- 
bicoherence for a different series of discharges in the con- 
figuration with i(a)/27T = 1.553. The elaboration of this 
figure is equivalent to that of Fig. [2l The mean line av- 
erage density at the transition is about 50% lower than 
with the previous configuration. Another difference with 
the previous result is that in these discharges, the tran- 
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FIG. 2: Series of 25 discharges in the configuration with 
i(a)/27r — 1.630. Total auto-bicoherence of complex Doppler 
reflectometry data, (a) Density profiles in the L and H phases, 
(b) Position p and k± of the measurement channels in the L 
and H phases, (c) Mean auto-bicoherence vs. time and nor- 
malized radius, p. (d) Mean auto-bicoherence vs. time and 
k±. (e) Mean auto-bicoherence vs. p and k± (averaged over 
— 10 < At < 10 ms). Radii and k± values of Figs, (c) - (e) 
correspond to the H phase. The vertical dashed line indicates 
the L-H transition; the horizontal dashed lines indicate the 
approximate position of the E r shear layer in the H phase. 
The statistical error level is about 0.006. 



sition at At — is not into the H-mode, but rather into 
the intermediate (I) phase between the L-mode (At < 0) 
and the H-mode (starting much later). We conjecture 
that this situation is due to a weaker zonal flow, inca- 
pable of producing the final transition into the H phase. 
This would be consistent with the observation of the 
predator-prey relation between turbulence and flows, as- 
sociated with wave numbers in the range k± ~ 6 — 12 
cm -1 , as confirmed by spectral analysis [7j. 

The location of the auto-bicoherence peak appears to 
be somewhat further inward than in the previous section 
< 0.76). 

Discussion and conclusions - In this work, we have 
used Doppler reflectometry measurements to show how 
the bicoherence evolves across the L-H transition at TJ- 
II. It is the first time that turbulent fluctuation am- 
plitudes and perpendicular (zonal) flow oscillations and 
their non-linear interactions have been observed simulta- 
neously in the region of the plasma where the transport 
barrier of the L-H transition forms. The observation 
of the bicoherence does not appear to depend very sen- 
sitively on the magnetic configuration. A surprisingly 
clear picture emerges that is consistent with established 
ideas about the interaction between zonal flows and tur- 
bulence [22| in the important formation phase of the H- 
mode transport barrier. 

In the case of the fast L-H transition, we have shown 
that the auto-bicoherence of the Doppler reflectometer 
signal is significant in a time window with a length of 
several tens of ms around the L-H transition time [l8| . 
at a particular radial position (p < 0.8, inward from the 
radial electric field shear layer in H mode), and that a 
specific range of perpendicular wave numbers is involved 
(k ± ~8- 11 cm" 1 ). Here, p s = ^2m~T e /eB = 0.1 - 0.2 
cm, so that this range corresponds to k±p s =0.8 — 2. 

A similar analysis in another set of discharges corre- 
sponding to a different magnetic configuration, in which 
the transition was characterized by an intermediate phase 
separating the L and H phases, the bicoherence appeared 
only after the L-I transition was made, and lasted several 
tens of ms afterward into the I phase. The perpendicular 
wave numbers involved (k± ~ 6 — 12 cm -1 ) were very 
similar to the previous case, while the radial location 
of the bicoherence peak was somewhat further inward 
(p<0.76). 

In both cases, the auto-bicoherence was due to high 
frequencies (of the order of 1-2 MHz) interacting with 
rather low frequencies (< 50 kHz), as expected for a zonal 
flow [13]. With the available data, it is not feasible to 
clearly resolve the low frequency component. 

Putting these results into context, we note that this is 
the first time wave number resolved bicoherence during 
L-H transitions is reported. This important result indi- 
cates which spatial scales are important in the zonal flow 
coupling producing the L-H transition. 
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FIG. 3: Series of 26 discharges in the configuration with 
l(o)/2ti = 1.553. Total auto-bicoherence of complex Doppler 
reflectometry data, (a) Density profiles in the L and I phases, 
(b) Position p and k± of the measurement channels in the L 
and I phases, (c) Mean auto-bicoherence vs. time and nor- 
malized radius, p. (d) Mean auto-bicoherence vs. time and 
k±. (e) Mean auto-bicoherence vs. p and k± (averaged over 
< At < 20 ms). Radii and k± values of Figs, (c) - (e) 
correspond to the I phase. The vertical dashed line indicates 
the L-I transition; the horizontal dashed line indicates the 
approximate position of the E r shear layer in the I phase. 
The statistical error level is about 0.006. 
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